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ABSTRACT 

Context. The two Oe stars HD 45 314 and HD 60 848 have recently been found to exhibit very different X-ray properties: whilst 
HD 60 848 has an X-ray spectrum and the emission level typical of most OB stars, HD 45 314 features a much harder and brighter 
X-ray emission, making it a so-called y Cas analogue. 

Aims. Monitoring the optical spectra could provide hints towards the origin of these very different behaviours. 

Methods. We analyse a large set of spectroscopic observations of HD 45 314 and HD 60 848, extending over 20 years. We further 
attempt to fit the Hff line profiles of both stars with a simple model of emission line formation in a Keplerian disk. 

Results. Strong variations in the strengths of the Hff, H/?, and Hei 4 5876 emission lines are observed for both stars. In the case 
of HD 60 848, we find a time lag between the variations in the equivalent widths of these lines, which is currently not understood. 
The emission lines are double peaked with nearly identical strengths of the violet and red peaks. The Ho- profile of this star can be 
successfully reproduced by our model of a disk seen under an inclination of 30°. In the case of HD 45 314, the emission lines are 
highly asymmetric and display strong line profile variations. We find a major change in behaviour between the 2002 outburst and the 
one observed in 2013. This concerns both the relationship between the equivalent widths of the various lines and their morphologies 
at maximum strength (double-peaked in 2002 versus single-peaked in 2013). Our simple disk model fails to reproduce the observed 
Hu line profiles of HD 45 314. 

Conclusions. Our results further support the interpretation that Oe stars do have decretion disks similar to those of Be stars. Whilst 
the emission lines of HD 60 848 are explained well by a disk with a Keplerian velocity field, the disk of HD 45 314 seems to have a 
significantly more complex velocity field that could be another signature of the phenomenon that produces its peculiar X-ray emission. 

Key words. Stars: early-type - stars: emission-line. Be - stars: individual: HD 45 314 - stars: individual: HD 60 848 - stars: winds, 
outflows 


1. Introduction 

By definition. Be stars are non-supergiant B stars that have at 
least once displayed Balmer line emission in their spectrum. Be 
stars are usually considered to be very rapidly rotating and non- 
radially pulsating B stars surrounded by an outwardly diffusing 
Keplerian disl£] About one fifth of the Galactic B stars belong to 
this category, but the strongest incidence of the Be phenomenon 
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1 In addition to the ‘e’ qualifier, several variants such as ‘(e)’, ‘pe’, 
‘[e]’ or ‘e+’ have been introduced. See Table 3 of Sota et al. J201 lb for 
a detailed definition of these variants. 


is found in the B1-2 subtypes (for recent reviews on Be stars, see 
Neiner l201 T, Rivinius et al. 12013 ). 

Conti & Leep d 19741 1 identified a small group of O-type stars 
- the so-called Oe stars - that seem to form an extension of the Be 
phenomenon towards higher stellar temperatures. Indeed, these 
stars display emission lines of the H i Balmer series, as well as of 
other elements, such as He i and Fe n, but do not exhibit conven¬ 
tional Of emission lines, such as He n A 4686 and N hi A 4634-40. 
Negueruela et al. (120041 note that previous spectral classifica¬ 
tions of Oe stars were probably too early because of infilling 
of the photospheric He i absorptions by circumstellar emission. 
These authors thus argued that most Oe stars have true spec¬ 
tral types in the range 09 - B0, with the notable exceptions 
of the 07.5 Die star HD 155 806 (Negueruela et al. 120041 and 
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HD 39 680, which was assigned an 06 V:[n]pe var spectral type 
by Sota et al. (1201 Ik 

Whilst the disks of some Be stars have been resolved in op¬ 
tical long baseline interferometry (see Rivinius et al. 120131 and 
references therein), the more distant and hence fainter Oe stars 
are beyond the capabilities of current interferometry facilities. 
Alternatively, the presence of a disk-like outflow can also be in¬ 
ferred from the depolarization of the continuum polarization due 
to line emission. The lack of such a signature in some Oe stars 
led Vink et al. (120091 to question the existence of disk-like struc¬ 
tures around these objects. Indeed, their linear spectropolarimet- 
ric observations revealed a clear depolarization effect across the 
Her emission line only in one case out of four Oe stars with Her 
emissior0 


The spectra of Be stars are strongly variable over a broad 
range of timescales. For instance, the relative strength of the 
peaks of double-peaked emission lines varies in so-called V/R 
(violet-to-red) irregular cycles that can last between a few weeks 
to decades with a mean around seven years (Neiner 12011 , 
Rivinius et al. 120131 ). The V/R cycles are usually attributed to 
precessing one-armed density waves in the disk (Okazaki 199 1). 
The precession periods of these density waves are two orders of 
magnitude longer than the orbital periods of the particles in the 
disk. Even more spectacular, the disk as a whole may dissipate 
and form anew (e.g. Draper et al. 120141 . 


With the exception of our previous study of the targets of the 
present paper and some scarce reports based on a limited number 
of data (Rauw et al. 120071 Copeland & Heard fl963l Divan et al. 
1983 1 Gamen et al.1 2012 ). the variability of Oe stars has received 
little attention so far. In the present paper, we study the variabil¬ 
ity of the emission lines in the spectra of HD 45 314 (PZGem) 
and HD 60 848 (BN Gem) classified as B0 IVe and 09.5 IVe, re¬ 
spectively, by Negueruela et al. (120041 ]] on various time scales 
from a few days up to several decades. Vink et al. (l2009t> report a 
spectropolarimetric detection of a disk in the case of HD 45 314, 
but not for HD 60 848. 


Both stars were recently studied in X-rays with XMM- 
Newton, revealing very contrasting X-ray properties (Rauw et 
al. 12013b . On the one hand, HD 60848 was found to display an 
X-ray emission similar to that of normal O-type stars, in terms 
of its X-ray luminosity and the temperature of the hot plasma. 
On the other hand, HD 45 314 exhibits an unusually hard and 
bright X-ray emission that makes it a new member of the cate¬ 
gory of so-called y Cas analogues. The y Cas analogues are de¬ 
fined by their remarkable X-ray properties (Smith et al. [20121 ). 
which include a ratio jyy ~ 10 fi , i.e. a factor 10 higher than for 
normal OB stars but at least an order of magnitude lower than 
for Be high-mass X-ray binaries, variations on timescales from 
minutes to months, and an X-ray spectrum that is dominated by 
a very hot plasma with kT > lOkeV. The origin of the y Cas 
behaviour currently remains unknown, but the most promising 
scenarios include a wind-disk interaction (Smith et al. 120121) or 


2 It must be stressed, though, that some other objects with known 
large-scale structures in their winds (e.g. 9 1 Ori C) were not detected 
either (see also Naze & Rauw l20l0ll . 

3 In the work of Conti & Leep 11974) . HD 45 314 and HD 60 848 
were classified as 09?pe and 08 V?pe, respectively. More recently, 
Sota et al. ( 2011 ) have assigned spectral types 09: npe and 08: V: pe 
to HD 45 314 and HD 60 848, respectively. Since Conti & Leep’s *?’ 
and Sota et al.’s qualifiers both mean uncertainties in the spectral 
subtype, both papers actually give the same spectral type for HD 60 848 
and only differ in the one for HD 45 314 in the ‘n’ qualifier (which 
means that the lines are broadened). 


perhaps accretion from the Be decretion disk onto a compact 
companion (White et al. 1982) . 

The goal of the present study is to characterize the spectro¬ 
scopic variability of these two Oe stars on a variety of timescales. 
For this purpose, we have conducted several dedicated observing 
campaigns, but we also retrieved a large number of spectra from 
various archives. 


2. Observations 

The data used in the present work come from a variety of ob¬ 
servatories and equipments (see Table|T]for details). The largest 
part of our data set was obtained via a collaboration between am¬ 
ateur and professional astronomers. Indeed, Be stars are promi¬ 
nent targets for amateur astronomer spectroscopists, and several 
collaborations have previously taken place between professional 
and amateur astronomers (e.g. Neiner [2011 ). In our case, most 
of the amateur data were taken during a four-month campaign 
(Eversberg [20111 ) conducted at the 50 cm Mons telescope located 
at the Observatorio del Teide and operated by the Instituto de 
Astrofisica de Canarias (IAC) in Tenerife. A LHIRES III spec¬ 
trograph was installed on the Mons telescope and operated be¬ 
tween early December 2008 and the end of March 2009. 

The main goal of this campaign was to monitor the colliding 
wind binary WR 140 during its periastron passage (see Fahed 
et al. [2011 ). Since WR 140 was observable for most of the run 
only close to either sunset or sunrise, additional targets were ob¬ 
served, including HD 45 314 and HD 60 848 (Morel et al. lWD . 
The wavelength coverage was between 6360 and 6950 A with 
a reciprocal dispersion of ~ 0.34 A pixel The data were re¬ 
duced using the IRAF software. To correct for some drifts in 
the instrumental set-up, we re-aligned the wavelength-calibrated 
spectra by means of the diffuse interstellar band at 6613.910 A or 
by means of strong telluric lines (e.g. 6892.3 A). After rejecting 
some spectra that were affected by strong ghosts due to rema- 
nence of the CCD after Ne-lamp exposures, we were left with 
35 usable spectra of HD 45 314 and 40 spectra of HD 60 848. 
The HD 60 848 time series provides a rather uniform coverage 
over the full duration of the campaign, whilst there is a gap in the 
middle of the campaign for the HD 45 314 time series. A few ad¬ 
ditional amateur spectra of HD 60 848 were retrieved from the 
BeSS database (Neiner [2011 1, operated at LESIA, Observatoire 
de Meudon, France]]] 

The rest of the data used in this paper were taken with 
professional instruments. Dedicated campaigns were conducted 
between February 1997 and September 2010 with the Aurelie 
spectrograph (Gillet et al. [1994 1 at the 1.52 m telescope of the 
Observatoire de Haute Provence (OHP, France). Various set-ups 
(both in terms of wavelength coverage and resolving power) 
were used. Between May 1999 and March 2002 echelle spec¬ 
tra were taken with the FEROS instrument (Kaufer et al. 11997 ) 
mounted on the 1.5 m telescope at La Silla (ESO, Chile). The 
data were reduced using the MIDAS software along with private 
software tools designed for the reduction of spectroscopic data 
taken with these instruments. Part of the Aurelie and FEROS 
data (up until 2002) were already used in the study of Rauw 
et al. (12007) . Another dedicated campaign was conducted be¬ 
tween September 2013 and February 2014 on the TIGRE tele¬ 
scope (formerly known as the Hamburg Robotic Telescope, 
Hempelmann et al. 120051 Mittag et al. 120111 Schmitt et al. 
12014 ) installed at La Luz Observatory (Guanajuato, Mexico). 


4 http://basebe.obspm.fr 
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Table 1 . Journal of the observations used in this work. 


Instrument 

Epoch 

Resolving power 

Wavelength domain 

Number of spectra 

HD 45 314 HD60848 

Aurelie @ 1.5m OHP 

Feb. 1997 

20000 

6510 - 6710 A 

11 

6 


Nov. 1998 

30000 

6500 - 6620 A 

7 

6 


Nov. 1998 

30000 

4795 - 4925 A 

1 

1 


Sep. 2000 

10000 

4460-4900 A 

3 

2 


Sep. 2001 

10000 

6350-6770 A 

3 

- 


Sep. 2001 

10000 

4460-4900 A 

1 

- 


Jan. 2008 

10000 

6350-6770 A 

1 

1 


Jan. 2008 

10000 

5490-5920 A 

3 

3 


Sep. 2008 

10000 

4460-4900 A 

1 

- 


Dec. 2009 

20000 

4450-4670 A 

1 

5 


Dec. 2010 

20000 

4450-4670 A 

- 

1 

FEROS @ 1.5m ESO 

May 1999 

48000 

3900-7100 A 

10 

10 


May 2000 

48000 

3900-7100 A 

4 

6 


May 2001 

48000 

3900-7100 A 

3 

3 


Mar. 2002 

48000 

3900-7100 A 

3 

3 

FEROS @ 2.2m ESO 

Jan.-Feb. 2006 

48000 

3900-7100 A 

3 

- 

RGO@ AAT 

Dec. 1994 

3500 

3700-6000 A 

1 

1 

IDS @ INT 

Jan. 1998 

4500 

4620-5020 A 

- 

2 


Oct. 2001 

3500 

3950-4880 A 

1 

1 


Jan. 2002 

4500 

6100-7100 A 

1 

1 

UVES @ UT2 

Jan. 2001 

54000 

3280-4560 A 

- 

1 


Jan. 2001 

87000 

4720-6830 A 

- 

1 

Elodie @ 1,9m OHP 

Nov. 2004 

42000 

3905-6810 A 

1 

1 

LHIRES III @ Mons 

Dec. 2008 - Mar. 2009 

3500 

6360-6950 A 

35 

40 

Sophie @ 1.9m OHP 

Mar. 2012 

40000 

3870-6940 A 

1 

- 

Coralie @ Euler 

Mar. - April 2012 

55000 

3850-6890 A 

2 

1 

FIES @ NOT 

Jan. 2011 

46000 

3700-7300 A 

1 

2 


Dec. 2012 

46000 

3700-7300 A 

1 

1 


Jan. 2013 

46000 

3700-7300 A 

3 

- 

HEROS @ TIGRE 

Sep. 2013 - Feb. 2014 

20000 

3500-8800 A 

10 

9 


The TIGRE features the refurbished HEROS echelle spectro¬ 
graph (Kaufer [1998 ) and is operated in a fully robotic way. 
We also include the echelle spectra taken with the Coralie and 
Sophie spectrographs, respectively, at the Swiss 1.2 m Leonhard 
Euler Telescope (La Silla, Chile) and the 1.93 m telescope at 
OHR which we obtained in March and April 2012 to support 
our X-ray observations of the Oe stars (see Rauw et al. 120131 . 

To complement our data, we also extracted spectra from a 
number of archives. We retrieved echelle spectra from the Fibre- 
fed Echelle Spectrograph (FIES) at the 2.56 m Nordic Optical 
Telescope at the Observatorio del Roque de Los Muchachos 
(La Palma, Canary Islands, Spain). These data were reduced 
with the IRAF software. Long-slit spectra were retrieved from 
the archives of the Intermediate Dispersion Spectrograph (IDS) 
at the Cassegrain focus of the 2.5 m Isaac Newton Telescope 
at the Observatorio del Roque de Los Muchachos and from 
the Royal Greenwich Observatory spectrograph at the 3.9 m 
Anglo Australian Telescope at the Anglo Australian Observatory 
(Siding Spring, Australia). The long-slit spectra were reduced 
and calibrated using the MIDAS software. Finally, calibrated 
echelle spectra were taken from the archives of the Elodie in¬ 
strument at the 1.93 m OHP telescope (Moultaka et al. 120041 , 
as well as from the UVES spectrograph at the ESO Very Large 
Telescope (Cerro-Paranal, Chile). 

For those data sets where we had observations of fast- 
rotating standard stars at different airmasses, we used the latter 
to build a template file of the telluric absorption lines. In all other 


cases, we used the telluric tool within IRAF along with the 
list of telluric lines of Hinkle et al. (12000b to remove the telluric 
lines in the Hei A 5876 and Her regions. Given the strengths of 
the Her emissions we are dealing with, the impact of the telluric 
lines on the equivalent widths is rather modest, well below the 
uncertainties due to the normalization. The impact of the correc¬ 
tion is somewhat greater for the weaker He i A 5876 lines. 

Figure]!] illustrates some important wavelength domains in 
the spectrum of HD 60 848 at a time when the circumstellar 
emission was at a low level. Beside the Balmer and Hei cir¬ 
cumstellar emission lines, we note the presence of a number of 
photospheric absorption lines, associated with Hen, Siiv, and 
N hi, which are typically observed in late-type O-stars. Figure[2] 
provides the same information for HD 45 314. In addition to 
these lines, the spectrum of the latter star also displays a num¬ 
ber of Fe ii emission lines. The lack of a strong photospheric 
He ii 4 4542 absorption indeed indicates that the central star of 
HD 45 314 must have a later spectral type than HD 60848. 

To evaluate the projected rotational velocity of our targets, 
we have extracted the line profile of the Hen 44686 absorp¬ 
tion in the low-emission state (HEROS spectra from December 
2013 for HD 60 848 and FIES spectrum from January 2011 for 
HD 45 314). Indeed, for both stars. Hen 44686 is the strongest 
absorption line that is free of emission from the circumstellar 
disk and is relatively free of blends with other lines. We then ap¬ 
plied the Fourier transform method (Gray [20051 Simon-Dfaz & 
Herrero [2007 ) to derive v sin/'. For HD 60 848 and HD 45 314, 
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Fig. 1 . Important spectral regions of HD 60 848 as observed with the HEROS spectrograph in December 2013 when the global 
emission level of the star was low. 


we obtain (230 ± 5)kms 1 and (210 ± 10)kms , respec¬ 
tively. These values are in good agreement with those obtained 
by Penny (fl996b : 231 km s ! for HD 60848 and 243 km s" 1 
for HD 45 314, via a cross-correlation method of IUE high- 
resolution spectra. Our value of HD 45 314 is slightly more un¬ 
certain because the blue wing of the He n 4 4686 line is blended 
with On /14676. In both cases, the derived value of v sin i is 
likely a lower limit. Indeed, as a result of gravity darkening, 
the He n absorption most likely forms nearer to the hotter stellar 
poles and is less sensitive to the cooler equatorial regions that 
rotate faster. 

3. Measurements 

Because we are dealing with a rather heterogeneous data set, 
combining data from different instruments and covering differ¬ 
ent wavelength domains at different spectral resolutions, we first 
need to decide what quantities we wish to measure. For the 
purpose of the present analysis, we focus on the radial veloci¬ 
ties (RVs), equivalent widths (EWs), and violet over red (V/R) 
ratio^] of those spectral lines that are most often covered by 

5 We define the V/R ratio as the ratio between the continuum- 
subtracted line intensities of the violet and red emission peaks. 


our observations. For HD 60 848, we obtain RVs for the Hen 
AA 4542, 4686, and 5412 absorption line^J For HD 45 314, we 
measure the RVs of the He ii 44686 absorption, as well as of the 
highest emission peak of the Hu line. He ii lines are intrinsically 
broader than He i or metallic lines, for example, and are there¬ 
fore subject to larger RV uncertainties. However, in the present 
case, the He n lines are in fact the only relatively strong and rela¬ 
tively isolated absorption lines that are not affected by emission 
from circumstellar material, and are thus suitable for RV mea¬ 
surements. For both stars, V/R ratios and EWs are measured for 
the Hj3, Her, and He i AA 5876,6678 lines. The EWs are measured 
over fixed wavelength windows. For HD 45 314, we use 4850 - 
4875, 6530 - 6600, 5865 - 5885, and 6670 - 6689 A for the Hj3, 
Her and Hei 445876, 6678 lines, respectively. For HD 60 848, 
we adopt the same windows, except for the Her line where we 
instead use 6550 - 6575 A. The uncertainties on the EWs were 
estimated using the formalism of Vollmann & Eversberg ( 20061) . 


6 We adopted the rest wavelengths 4541.59, 4685.68, and 5411.52 A 
from Underhill d 1994b . 
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Fig. 2. Same as Fig. jTJbut for HD 45 314 as observed with the FIES spectrograph in January 2011 when the global emission level of 
the star was low. 


4. Results 

4.1. HD 60 848 

Figure[3] illustrates the long-term behaviour of the Ha, H/3, and 
He i 45876 emissions in the spectrum of HD 60 848. The first 
maximum of EW(Hcr) and the subsequent decay have already 
been reported in Rauw et al. (120071) . Figure[3] reveals another 
episode of maximum Her emission around HJD 2454500. By 
chance, the Mons campaign covers the decline of this episode 
and reveals a rather monotonic and fast drop of EW(Hcr) (see 
Fig. [4]). We have searched for periodicities in the EW(Hcr) vari¬ 
ations using both the Fourier-method for uneven sampling of 
Heck et al. dl985b . modified by Gosset et al. (120011 ). and the 
string-length method of Lafler & Kinman (i 1965b , the latter be¬ 
ing more appropriate for highly asymmetric variations. In this 
way, we find that the available observations suggest a recurrence 
timescale of about 2775 days. 

There is generally a good correspondence between the be¬ 
haviour of the EWs of Her and Hei 45876, although the varia¬ 
tions in the latter lag a bit behind those of the former. This lag 
is clearly seen in Fig. [5] where we show the EWs of both lines 
whenever they are available simultaneously. 


The situation is somewhat different for H/3. For those ob¬ 
servations where we have simultaneous access to the Her and 
H/3 lines, we have again plotted the equivalent width of the for¬ 
mer line versus that of the latter. The result is shown in Fig. [6] 
When the emission lines have their lowest level, the H/3 line is 
dominated by the photospheric absorption, leading to a positive 
equivalenth width, whilst Her is always seen in pure emission 
in our data. When the emission level increases, the equivalent 
widths of H/3 increase at a faster rate than EW(H<r) (see Fig. [6}. 
This trend continues up to EW(H/j) - -0.8 A, and EW(Hcr) 
^ -6.5 A. Beyond that point, the slope of the relation levels off 
and EW(Ho') now increases much faster than EW(H/j). The de¬ 
cline of the emission level occurs very rapidly. 

The radial velocities of the Hen absorption lines display 
some dispersion that might a priori be due to variability on rela¬ 
tively short timescales. For the Hen 444542, 4686, and 5412 
lines, we obtain mean RVs of 35.4 + 11.4, 24.0 + 5.7, and 
25.7 ± 7.5 km s 1 , respectively. Here the quoted dispersions cor¬ 
respond to standard deviations about the mean values. We have 
searched for periodicities using the methods introduced above, 
but no consistent periodicity was found. Boyajian et al. (120071 ) 
used a bisector method to measure the velocities of the emis¬ 
sion lines of HD 60 848 during an observing run in December 
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Fig. 3. Long-term variations in the equivalent widths (in A) of 
the H<r, H/j, and He i A 5876 lines in the spectrum of HD 60 848. 
The lower panel yields the meaning of the different symbols. For 
each of the three upper panels, the lower horizontal axis yields 
the date in HJD -2450000, whilst the upper horizontal axis 
yields the epoch expressed as Gregorian calendar years. Typical 
uncertainties on the EWs are 0.25, 0.25, and 0.20 A for the Her, 
H/i, and Her 45876 lines, respectively. The Coralie spectrum 
was taken simultaneously with the XMM-Newton observation 
discussed by Rauw et al. (120131) . 
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Fig. 5. EW(He i A 5876) as a function of EWfHcr) for HD 60 848. 
The dashed line connects the data points taken during the high 
emission cycle between 1999 and 2004. The symbols have the 
same meaning as in Fig. [3] 
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Fig. 6. Same as Fig. [ 5 ] but for EW(HyS) as a function of EW(Hcr). 


Fig. 4. Zoom on the trends observed during the intense moni¬ 
toring of the Mons campaign for the equivalent widths and V/R 
ratio of the Her line in the spectrum of HD 60 848. Typical un¬ 
certainties are 0.25 A for the EW and 0.01 for the V/R ratio. The 
insets in the top panel illustrate the first (18 December 2009, left 
inset) and last (14 March 2010, right inset) normalized spectra 
around the Her line as observed during the Mons campaign. 


2000. They found possibly recurrent RV variations on timescales 
of three to four hours. These timescales are too short to be as¬ 
sociated with orbital motion in a binary system, but could be 
due to non-radial pulsations, as suggested by McSwain et al. 
(120071 ). Whilst our sampling is not suited to checking for the 
existence and stability of the periodicities reported by Boyajian 
et al. ( 120071 ). we note that the absorption lines indeed display line 
profile variability that might be at the origin of the RV variations. 
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The V/R ratio of the Her emission remains very close to 
unity. The standard deviation about the mean value (1.00) of 
our 107 measurements amounts to 0.06. Our search for peri¬ 
odicities yields no significant result. The V/R ratios of H/3 and 
He i /l 5876 behave similarly, although in these cases, the am¬ 
plitude of variation increases when the overall level of emission 
goes down (which is as expected in view of the definition of 
V/R). The separation between the emission peaks of the Her line 
changes from about 210 - 220 km s 1 at minimum emission to 
150 - 160l<ms 1 at maximum emission level. In a Keplerian 
disk, where the separation between the peaks is interpreted as 
the value of 2 VK ep sin i at the outer rim of the disk (Horne & 
Marsh [19861) . such a variation indicates that the outer radius of 
the disk varied by a factor ~ 1.9 between minimum and maxi¬ 
mum emission. We come back to this point in Sect. 0 

4.2. HD 45314 

Figure[7]shows the long-term behaviour of the equivalent widths 
of the Ha, H/3, and Hei A 5876 emissions in the spectrum of 
HD 45 314. These lines are always seen in pure emission in 
our data, even at those times of minimum global emission level. 
Copeland & Heard ( 119631 ) reported the absence of H/3 and Hy 
emission in 1962. Our data contain no such episode, at least not 
for H/3. 


Epoch 
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Fig. 7. S ame as Fig. [ 3 ] but for HD 45 314. Typical uncertainties 
on the EWs are 0.50, 0.25, and 0.20 A for the Ha, H/3, and He 1 
A 5876 lines respectively. We have also included the EW(Ha) 
measurement of Silaj et al. (12010b . The second Coralie spec¬ 
trum was taken simultaneously with the XMM-Newton observa¬ 
tion discussed by Rauw et al. (12013b . 


Figure[7] reveals up to three episodes of high emission, al¬ 
though with different behaviours in the three lines. During the 
first maximum (HJD ~ 2452500), all three lines are in strong 
emission. The Mons campaign seems to cover another episode 
of slowly increasing Ha emission (see Fig. [8}, roughly 3200 days 


after the first event. Unfortunately, we lack contemporaneous 
data in the other two lines. In our data a clear change in be¬ 
haviour occurs for the third episode of strong emission af¬ 
ter HJD 2456000. Whilst the H/3 and Hei 45876 lines both 
strongly increase, reaching their level of the first maximum, 
Ha does not follow the trend and remains at a rather moderate 
level of cmissiorQ The behaviour of HD 45 314 thus appears to 
change from one outburst to another. 



4600 4820 4840 4860 4880 4900 4920 



HJD - 2 450 000 


Fig. 8. Same as Fig. [ 4 ] but for HD 45 314. Typical uncertainties 
are 0.50 A for EW and 0.01 for the V/R ratio. The insets in the 
lower panel illustrate the first (17 December 2009, left inset) and 
last (22 March 2010, right inset) normalized spectra around the 
Ha line as observed during the Mons campaign. 


The general behaviour of HD 45 314 during the 2013 event 
indeed seems different from the 2002 spectroscopic outburst. 
This is illustrated in Fig. [9] In the 2002 event, the emission 
lines displayed similar morphologies at minimum and maxi¬ 
mum, whilst these morphologies drastically changed between 
minimum and maximum in 2013. In December 2013, i.e. near 
maximum, all three lines were single-peaked. The Ha emission 
remained single-peaked at least until February 2014. Whilst the 
maximum Ha intensity reached roughly the same value as in the 
previous spectroscopic outburst, the big difference is the width 
of the line. In fact, in December 2013, the Ha line at its maxi¬ 
mum intensity is significantly narrower than previously. 

The V/R ratios of the emission lines in the spectrum of 
HD 45 314 undergo large and rather slow variations between 
0.82 and 1.56 for Ha, between 0.60 and 2.07 for H/3, and be¬ 
tween 0.66 and 1.88 for He 1 A 5876 (see Fig. 10 1 . For those spec¬ 
tra where a single peak was observed, the measurements were 
discarded. Since the morphology of the Ha line does not always 
allow the strengths of the peaks to be easily measured, we have 
also used the RV of the dominant peak as a proxy of the V/R 
variations of this line. A Fourier analysis of the time series of 
V/R measurements yields low-significance peaks at medium and 
long timescales. For the H/3 line, the V/R variations are domi¬ 
nated by a timescale of ~ 160 days. For the RVs of the Ha peak 
and the V/R of Hei 45876, we find a dominant timescale of 


7 We have checked that none of the Ha measurements used here is 
affected by saturation of the CCD. Some of the FIES @ NOT data suffer 
from saturation in this line, but they were discarded from our analysis. 
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Fig. 9. Line profile variations of HD 45 314 between minimum 
and maximum emission at different epochs. The top row illus¬ 
trates the variations in the spectrum between 7 May 2000 (black 
solid line, HJD 2451 672.5) and 5 March 2002 (red dotted line, 
HJD 2 452 338.5). The bottom row yields the spectra as observed 
on 28 September 2013 (black solid line, HJD 2456 564.0) and 
20 December 2013 (red dotted line, HJD 2 456 646.9). 

Epoch 


1995 3000 3005 3010 



HJD - 3 450 000 

Fig. 10. Variations in the V/R ratio of the Her, H (3, and Hei 
4 5876 emission lines in the spectrum of HD 45 314. The sym¬ 
bols have the same meaning as in Fig. [7] 


~ 420 days in both cases. Finally, the V/R ratio of Her is modu¬ 
lated on a timescale of 2500 days. The latter is also present in the 
periodograms of the V/R ratios of the H/3 and He i A 5876 lines, 
and of the RV of the dominant Her peak. We thus conclude that 


two possible timescales are 420 and 2500 days, although their 
detections are not statistically significant. 

Unlike what we found for HD 60 848, in the case of 
HD 45 314, we observe no evidence for a systematic time delay 
between the variations in the EWs of Her and Hei 45876 (see 
Fig-Oj}: the EWs of both lines seem well correlated during the 
outbursts though with a shift between the 2002 and 2013 events. 
In contrast, the relationship between EW(Hcr) and EW(H/J) is 
more complex (see Fig.fTT}. At first sight the behaviour during 
the 2002 spectroscopic outburst is similar to what was seen in 
HD 60 848 (compare with Fig.|6|. However, as for the plot of 
EW(Hei 45876) versus EW(Hcr), the 2013 event is clearly off¬ 
set with respect to this behaviour. This is due to the significantly 
lower level of Her emission in 2013. 
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Fig. 11. EW(Hei 45876) as a function of EW(Hcr) for 
HD 45 314. The dashed line connects the data points taken 
around the 2002 spectroscopic outburst, whilst the dotted line 
connects the data from the 2013 outburst. The symbols have the 
same meaning as in Fig. [7] 


The RVs of the Hen 44686 absorption line yield a mean 
value of -1.8+14.3 kms _I . As for HD 60 848, we have searched 
for periodicities using the Fourier and string-length methods in¬ 
troduced above. No significant periodicity is found in our RV 
time series. Boyajian et al. (20071 ) found a significant RV shift 
between their data of HD 45 314 taken in October and December 
2000. As for HD 60 848, the RVs of Boyajian et al. (l2007i > refer 
to bisector velocities of the emission lines, and the shift reported 
by these authors could thus actually be due to changes in the 
shape of the emission line profiles. 


5. Discussion 

5.1 . Fitting the Ha profiles 

Assuming that both Oe stars feature a Keplerian decretion disk, 
we can attempt to quantify the properties of these disks at dif¬ 
ferent epochs. For this purpose, we designed a simple model 
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Fig. 12. Same as Fig.|TT] but for EW(H/3) as a function of 
EW(Hor). The dashed line connects the data taken around the 
2002 outburst, whilst the dotted line corresponds to the 2013 - 
2014 event. 
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that allows the emission profile from an optically-thick disk to 
be computed. Our numerical code is based on the approach of 
Hummel & Vrancken (120001) . which in turn follows the formal¬ 
ism of Horne & Marsh (11986b for computing accretion disks in 
cataclysmic variables. A very similar approach was adopted by 
Grundstrom & Gies (12006b and Catanzaro d2013b to compute 
synthetic Be emission line profiles and compare them with ac¬ 
tual observations. 

The disk is assumed to be in pure Keplerian rotation and to 
extend from the stellar surface out to an outer radius Rdisk- The 
particle number density in the disk is described as a power law 
of index a in radius (r) and a Gaussian profile in elevation (z) 
above the disk plane: 


n(r, z.) = n 0 



1 


2 \H(r) 


( 1 ) 


where H(r ) is the local scale height (see equations 2 and 4 of 
Hummel & Vrancken [20001 . At each position in the disk, the 
line-of-sight optical depth (r) is computed accounting for the 
line profile width due to the thermal velocity and the shear broad¬ 
ening (see equations 7, 8, and 9 of Hummel & Vrancken 12000b . 
The specific intensity along the line of sight is computed accord¬ 
ing to 


h = fii(Tdisk) [1 - exp(-r)] + I\ exp (-r), (2) 

where we assume an isothermal disk and adopt the Planck func¬ 
tion at Tdisk = 0.6 7/ii as the source function. Occultation of the 
stellar surface by the disk is accounted for by the 7* exp (—t) 
term, and occultation of the disk by the stellar body is also taken 
into account. Since we assume an axisymmetric disk, the syn¬ 
thetic line profiles are by definition symmetric about the central 
wavelength. We note that electron scattering is not accounted for 
in our model. 


For our computations, we have adopted 7/n = 32500 K, 
R* — 8 R 0 , and M* = 20 M Q which are the mean values of the 
stellar parameters of HD 45 314 and HD 60 848 as quoted by 
Vink et al. (12009b . With these numbers fixed, there remain four 
free parameters in the model: the orbital inclination /, the outer 
disk radius Rdisk- the power-law index a, and the particle number 
density at the inner edge of the disk no- 

To account for the underlying photospheric absorption line, 
we adopted the effective temperatures (31500 and 33400 K for 
HD 45 314 and HD 60 848, respectively) and surface gravity 
(logg = 3.92 for both stars) as quoted by Vink et al. (12009b . 
We then interpolated in the Ostar2002 grid of non-LTE plane- 
parallel TLUSTY models (Lanz & Hubeny 20031 for solar 
metallicity. The resulting spectrum was broadened with the pro¬ 
jected rotational velocities inferred above to obtain the photo- 
spheric spectra. 

In our comparison with the actual observations, we focused 
on the Her line. First we computed a grid of synthetic line profiles 
over a wide range of parameters (i from 15 to 75° with a step of 
5°, Rdisk between 2 and 30 R* by stpdf of 1 R„, a between 1.5 and 
3.5 with a step of 0.5, and no varying between 16 and 626cm 3 
with a step of 16 cm -3 ). For each star, we considered the ob¬ 
served Her profile at several epochs: 1999, 2000, and 2002. For 
each of these epochs, we averaged all the FEROS @ 1.5m ESO 
data that are available. We also included the average HEROS 
@ TIGRE spectrum of HD 60848 as observed in 2013 - 2014. 
For HD 45 314, we further considered the mean spectrum of the 
2011 FIES @ NOT dataset and the HEROS @ TIGRE spectrum 
taken on 20 December 2013 (see Fig.[9]t. 

HD 60 848 presents only small-amplitude V/R variations 
and thus displays quite symmetric line profiles. These proper¬ 
ties make this star well suited for application to our model. 
Fitting our model to the H<r line profiles from all four epochs 
indeed yields good results. As an example. Fig. [13] displays the 
fits for two epochs that roughly correspond to the most extreme 
values of EW(Hcr). The best-fit model parameters are listed in 
Table[2] As can be seen, all epochs yield an inclination of the 
disk near 30°. The inclination inferred from the 1999 data is 
slightly larger than at the other epochs, but this may not be sig¬ 
nificant. The best-fit inclination yields an equatorial rotational 
velocity of 460kms -1 , i.e. about 67% of the critical velocity 
for the stellar parameters adopted here. However, as pointed out 
above, this estimate could be biased by gravity darkening effects, 
because we have derived the v sin i value from the He 11 A 4686 
line, which is likely to arise from the polar regions where the 
rotation is slower. 

All observations are consistent with a density-law exponent 
a - 3.0. The outer radius of the disk varies from about 12R* 
at minimum emission state to 28 R* at the intermediate emission 
state observed in May 1999. Surprisingly, the spectrum with the 
strongest Htr emission (May 2000) does not correspond to the 
largest disk radius, but rather corresponds to the highest density 
no- We repeated the fitting by fixing the inclination to 30° and 
a to 3.0. The main changes with respect to the previous results 
are a reduction of the disk radius for the May 1999 epoch to 
(21.5 ± 1.7)R, associated with a 17% reduction of uq. The disk 
radius in May 1999 is now comparable to the value inferred for 
the May 2000 epoch (20.5 R,). Our results thus suggest that the 
high state in HD 60 848 corresponds to an increase in disk radius 
followed by an increase in the density in the disk. 

The best-fit parameters in Table[2 are broadly consistent with 
the properties inferred for classical Be stars. Indeed, « () values 
inferred by Hummel & Vrancken (120001 ) span a wide range from 
33 to 1005 cm -3 , and our results are well inside this range. The 
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Fig. 13. Fits of the Her line of HD 60 848 as observed in May 
2000 (left) and the winter season 2013-2014 (right) with our 
disk model. The model is shown by the dashed line. The cor¬ 
responding model parameters are given in Table[2] The failure 
to reproduce the line wings is probably because our synthetic 
models neglect electron scattering. 


Table 2. Best-fit parameters of the Her line profile of HD 60 848. 
The quoted uncertainties correspond to the 90% confidence in¬ 
terval. 


Epoch 

ni 

Rdisk (R*) 

n 0 (cm" j ) 

a 

May 1999 

36.7 ± 3.3 

27.9 ± 1.6 

420 ± 117 

2.9 ±0.2 

May 2000 

30.0 ± 3.0 

20.4 ± 1.1 

585 ± 31 

3.0 ±0.2 

Mar. 2002 

28.2 ± 3.4 

14.4 ± 2.7 

270 ± 94 

3.1 ±0.3 

2013-14 

29.2 ± 3.4 

12.6 ± 2.3 

199 ±92 

2.8 ± 0.4 


density-law exponents of Be disks are close to 3 (with a range 
from 1.6 to 3.6, Hummel & Vrancken 2000 ) as we also find 
for HD 60 848. The only parameter that stands out in the case 
of HD 60 848 is the outer radius of the disk. Be disks are usu¬ 
ally found to be smaller with most objects having R disk less than 
12 R* (Grundstrom & Gies 120061 Catanzaro 12013b and only a 
few objects seem to have disks as large as 18.5 R* (Hummel & 
V rancken 1 2000 ). 

The situation is more complex for HD 45 314. Indeed, the 
observed line profiles of this star are more asymmetric and dis¬ 
play significant V/R variations. Moreover, the line profile vari¬ 
ations not only concern the relative strengths of the peaks, but 
also result in bumps in the emission wings. Furthermore, the 
observed profiles display broad emission wings that suggest a 
velocity at the inner rim of the disk that exceeds the velocity at 

the stellar radius calculated from the adopted values of 

M* = 20 M 0 and R, = 8 Rq. All these facts suggest that the ve¬ 
locity field in the emission volume around HD 45 314 is more 
complex than the simple smooth Keplerian disk that we assume 
in our model. Unsurprisingly, our model fails to correctly re¬ 
produce the observations. The least disappointing results are ob¬ 
tained for the FIES@NOT spectra taken in 2011 (see Fig.[l4|. 
Still, this fit has much lower quality than what we achieved for 
HD 60 848. For the May 1999, May 2000, and January 2011 


observations, we obtain a ‘best-fit’inclination of 45° and a disk 
radius of - 22.5 R*. The index of the density law is found to 
be 2.5. However, this apparent consistency between the parame¬ 
ters does not hold for observations taken when the star is in high 
state. Indeed, the March 2002 data yield a ‘best-fit’ inclination 
of 60°, whilst the 2013 spectra suggest 25°. This situation essen¬ 
tially reflects our failure to achieve a decent fit of the observed 
profiles with our model. 




RV (km/s) 

Fig. 14. Fits of the Ha line of HD 45 314 as observed in March 
2002 (left) and January 2011 (right) with our disk model. The 
model is shown by the dashed line. 

The formal quality of the fit to the broad wing^jof the 2002 
spectra could be improved by artificially increasing the M,jR, 
ratio by a factor ~ 4. Such a high ratio is, however, at odds 
with the typical parameters of Population I O-type stars, and the 
agreement would not be any better if assuming typical masses 
and radii of subdwarf B-stars (Hebcr j 1986) . 

As pointed out above, our failure to correctly reproduce the 
line profiles probably points towards a significantly more com¬ 
plex velocity field than assumed in our model. Whilst the ob¬ 
served asymmetries of the line profiles can be explained, at least 
to first order, assuming a sinusoidal modulation of the disk den¬ 
sity, the changes in line widths and in the inferred disk inclina¬ 
tion are more challenging. Genuine inclination changes could 
stem from the torque from a neutron star companion if the 
disk were misaligned with the neutron star’s orbit. Martin et 
al. (1201II) show that in such a configuration, the disk could be 
truncated inside the neutron star’s orbit, preventing Type II X- 
ray outbursts. Such a disk would be warped and would precess 
on a timescale that is much longer than the neutron star’s orbit. 
However, we currently lack clear evidence of such a neutron star 
companion in HD 45 314 because there is no convincing signa¬ 
ture of RV variations due to binarity, and the X-ray properties of 
this star differ from those of classical high-mass X-ray binaries. 
Moreover, given the poor quality of some of the line profile fits, 
the evidence of inclination changes is rather weak. 

8 These broad wings are not artefacts due to normalization errors of 
the echelle spectra. In fact, the same features are seen in the Aurelie 
data, which are much easier to normalize. 
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Finally, we have tested whether a non-Keplerian velocity 
law could improve the situation. For this purpose, we com¬ 
puted models with a velocity law v rot (r) = (y) ^ 

this relation [3 = 0.5 corresponds to a Keplerian velocity field, 
whilst P = 1 corresponds to conservation of angular momentum. 
Neither (3 — 0 nor (3 — 1 allow improvement in the quality of the 
fit of the observed profiles. 

5.2. The time dependence of the spectra 

In the case of the B2IV-Ve star cj CMa, Stefl et al. (120031) report 
a time delay of several months between the increase in the con¬ 
tinuum flux and higher order Balmer lines, on the one hand, and 
the increase in Her, on the other hand. These authors note further 
that the path from quiescence to high-emission state is distinct 
from the return path in the Balmer decrement diagrams (Hy/H/J 
versus HcrjHfi). As a result, the star describes a hysteresis-like 
path that Stefl et al. (120031) interpret in terms of the build-up and 
decay of the disk, both of which progress radially outwards. In 
this scenario, the mass-weighted mean radius of the disk would 
be smaller during build-up than during decay. 

The time dependence of the equivalent widths of the Ha 
and H(3 emission lines in Fig. [6] indicates that in the case of 
HD 60 848, Her also seems to lag behind H/i, at least during the 
early stages of the outburst. At first sight this is reminiscent of 
the situation of co CMa. That the increase of the H/j line levels off 
whilst Her continues to increase could indicate that the excitation 
of the n = 3 level of hydrogen in the disk continues to progress 
outwards, once the region of the disk where the n — 4 level is 
excited has reached its maximum extension. Alternatively, when 
the disk becomes very dense, parts of it could become cooler 
leading to a reduction in the relative population of the upper lev¬ 
els of hydrogen and thus a slower increase in H(3. However, these 
scenarios are contradicted by Fig. [5] which indicates that Hei 
A 5876 lags behind the variations in Ha: the EW of He i A 5876 
continues increasing once EW(Ha) has reached its maximum, 
even though the excitation of the upper level of Hei A 5876 re¬ 
quires almost twice as much energy as the excitation of the n — 3 
level of hydrogen. 

The situation is again different for HD 45 314. No hysteresis¬ 
like loops were found, but the star changes its behaviour from 
one epoch to the next. Whether this is related to the status of 
HD 45 314 as a y Cas analogue (Rauw et al. 1201 31 is not clear. 
At least the transition between double- and single-peaked emis¬ 
sion lines is not restricted to y Cas analogues because such tran¬ 
sitions have been reported for other Be stars (e.g. Hanuschick 
1 986 ). To our knowledge, apart from y Cas (B0.5 IVpe) itself 
(Miroshnichenko et al. [2002 , Pollmann et al. 120141 . none of the 
currently known y Cas analogues have been monitored spectro¬ 
scopically on long timescales. In their long-term study of y Cas, 
Miroshnichenko et al. (120021 reported V/R variations with cy¬ 
cles of progressively increasing length (~ 4, 5, 7, and 10 years). 
The V/R ratios of the Ha and non-hydrogen lines in y Cas 
were found to vary non-simultaneously, with the Ha cycle be¬ 
ing longer. In HD 45 314, we do not find evidence of any such 
effect, since the Ha, H/l, and He i A 5876 lines all display similar 
timescales in their V/R variations. 

6. Summary 

In this paper we have reported analysis of an extensive set of op¬ 
tical spectra of the two Oe stars, HD 45 314 and HD 60 848. The 
morphology of the Ha emission of HD 60 848 is explained well 


by a Keplerian disk around the star, despite previous reports of 
a lack of a spectropolarimetric signature for such a disk (Vink et 
al. 2009 ). The emission lines of HD 60 848 undergo strong vari¬ 
ations in their strength, with H/l increasing first, followed by Ha, 
and finally by He i A 5876. This situation cannot be explained by 
simple considerations of the progressive excitation of the disk, 
and the origin of the time lags is currently unclear. 

Whilst HD 45 314 is the only Oe star in the sample of Vink et 
al. (120091) where a clear depolarization effect across the Ha line 
was found in the linear spectropolarimetry, our simple Keplerian 
disk model fails to reproduce the observed line profiles of this 
star. This star furthermore exhibits strong variations in the mor¬ 
phology and strength of its emission lines. Moreover, we found a 
change in behaviour between previous high-emission states and 
the most recent event (2013 - 2014). Currently, the reasons for 
this strange behaviour and the connections with the fact that this 
star is a y Cas analogue are unclear. Our X-ray observation of 
HD 45 314 was obtained when the star was in a low state. If the 
y Cas status is directly related to phenomena in the disk, as sug¬ 
gested by Smith et al. (120121) . one would expect the X-ray emis¬ 
sion to be different during a high state. We will thus continue 
our optical monitoring of this star with the goal of triggering an 
X-ray observation during its next outburst. 

In the sample of Vink et al. ( 120091 ), four Oe stars display a 
strong Ha emission (HD 45 314, HD 60848, HD 120678, and 
HD 155 806), but only HD 45 314 displays a strong depolariza¬ 
tion effect. One can thus speculate whether the positive detection 
of HD 45 314 might stem from one of its peculiarities: the strong 
V/R variations, hinting at a precessing density wave, or its sta¬ 
tus as a y Cas analogue. The first possibility is unlikely because 
the spectra of HD 120678 (08III:nep) collected by Gamen et 
al. (120121) over five years also reveal strong V/R variation^] 
for this star, although it shows no strong depolarization effect. 
Concerning the second possibility, we note that whilst nothing 
is known about the X-ray emission of HD 120678, the X-ray 
emission of HD 155 806 and HD 60848 are typical of normal 
OB stars (Naze et al. I20l0l Rauw et al. [20l3T ) and HD 45 314 is 
thus currently the only known y Cas analogue among the Vink 
et al. (i2009t sample. Therefore, we cannot rule out the possibil¬ 
ity that the strong depolarization effect in HD 45 314 might be 
related to its status as a y Cas analogue rather than to the mere 
presence of a disk. 
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